Magnetization-dependent T c shift in F/S/F trilayers with a strong ferromagnet 
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We have measured the superconducting transition temperature T c of Ni/Nb/Ni trilayers when the 
magnetizations of the two outer Ni layers are parallel (P) and antiparallel (AP). The largest difference 
in T c occurs when the Nb thickness is just above the critical thickness at which superconductivity 
disappears completely. We have observed a difference in T c between the P and AP states as large as 
41 mK - a significant increase over earlier results in samples with higher T c and with a CuNi alloy 
in place of the Ni. Our result also demonstrates that strong elemental ferromagnets are promising 
candidates for future investigations of ferromagnet/superconductor heterostructures. 

PACS numbers: 74.45.+C, 85.75.-d, 85.25.-j, 73.43. Qt 
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Heterostructures composed of ferromagnetic (F) and 
superconducting (S) materials have attracted much theo- 
retical and experimental attention due to the rich physics 
produced by the interplay between competing symme- 
tries of the order parameters Q. In an S/F bilayer the 
exchange field of the ferromagnet modulates the super- 
conducting order parameter as it decays inside the fer- 
romagnet over a very short distance. Kontos et al. 
used tunneling spectroscopy to observe the damped os- 
cillations of the order parameter by measuring the den- 
sity of states (DOS) for different thickness ferromagnets. 
Ryazanov et al. [j| observed 7r-state Josephson coupling 
in an S/F/S trilayer first by varying temperature, then 
later by varying the thickness of the ferromagnet 0] . Ear- 
lier, several groups 0, IS na d observed oscillations in 
the critical temperature T c of S/F bilayers as a function 
of the ferromagnet thickness cIf- Under ideal conditions 
T c oscillations arise from interference between the trans- 
mitted superconducting wave function through the S/F 
interface and the wave reflected from the opposite surface 
of the ferromagnet, although in some cases alternative 
explanations have been proposed : 8j. In many experi- 
ments, weakly ferromagnetic alloys were used in order to 
reduce the size of the exchange splitting in the conduc- 
tion band, E ex , and thus increase the penetration length 
£f for Cooper pairs, where £f = HvF/2E ex in the clean 
limit and vf is the Fermi velocity of the ferromagnet 

An alternative way to probe the influence of a ferro- 
magnet on a superconductor is to look for T c variations 
in an F/S/F trilayer structure based on the mutual ori- 
entation of the two ferromagnet magnetizations [ToL fill . 
This effect was observed [l2J and later reproduced |l3| 
in a Cui-zNi^/Nb/Cui-^Nij, system, where a weak fer- 
romagnet was used because it is "less devastating to su- 
perconductivity." The largest difference in T c observed 
between the antiparallel (AP) and parallel (P) states of 
the F-layer mutual magnetizations was only 6 mK when 
T c was 2.8 K. Unlike other experiments HHSHIE0 
that require the ferromagnet thickness to be comparable 
to however, a positive feature of this experiment is 
that the T c difference is predicted to persist even for thick 



F layers 0, 0] . Thus it proves advantageous in study- 
ing systems with strong elemental ferromagnets, which 
have extremely short values of £p . 

Experimental studies of F/S systems with strong fer- 
romagnets are of interest because they provide new chal- 
lenges to theory, which does not yet address the full com- 
plexity of the ferromagnetic state with its different DOS 
and vp of the majority and minority spin bands. Fur- 
thermore, pure elemental ferromagnets are in the clean 
limit, £p < If where If is the mean free path; this com- 
plicates use of the the popular Usadel equations normally 
applied to the dirty limit. We are motivated to work in 
the limit of thick ferromagnetic layers, in anticipation 
of future situations where superconducting order with 
spin-triplet symmetry is induced in a superconductor sur- 
rounded by ferromagnets with noncollinear magnetiza- 
tions [l^. (When dp 3> £f, the singlet component of 
the order parameter is completely damped.) Lastly, we 
wish to understand whether, in an F/S/F system with 
a strong ferromagnet, a large difference in T c between 
the P and AP states can be achieved, as envisioned in 
the proposals for a superconducting spin switch [Tfl lllj | . 
The Ni/Nb system has been shown to be a viable can- 
didate for experiments on F/S systems 0,0- I n this 
paper we will show that Ni/Nb/Ni trilayers exhibit a sig- 
nificant T c shift depending on the mutual orientation of 
the magnetizations of the two Ni layers. 

Sets of Ni(7) /Nb(4)/Ni(7) /Fe 50 Mn 50 (8) /Nb(2) multi- 
layers (all thicknesses are in nm) were directly deposited 
onto Si substrates by magnetically-enhanced triode dc 
sputtering in a high vacuum chamber with a base pres- 
sure in the low 10 -8 Torr and an Ar pressure of 2.0 ■ 1CP 3 
Torr. The Ni thickness of 7 nm was chosen to be much 
longer than £p, which we estimate to be 0.8 nm using 
2E ex = 0.23 eV and v F = 0.28 ■ 10 6 m/s for the major- 
ity band 0|- The purpose of the FeMn is to pin the 
magnetization direction of the top Ni layer by exchange 
bias [lij . The non-superconducting Nb capping layer 
protects the FeMn from oxidation. After deposition, the 
samples were heated to 180°C under vacuum, just above 
the blocking temperature of FeMn, and cooled in an ap- 
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FIG. 1: Critical temperature vs. Nb thickness for 
Ni(7)/Nb(d s )/Ni(7)/Fe 50 Mn 5 o(8)/Nb(2) samples (all thick- 
nesses are in nm). Different symbols represent different sput- 
tering runs. The solid line represents the theoretical fit. Inset: 
Schematic cross-section of the samples. 

plied field of 200 Oe in the plane of the multilayer. This 
procedure pins the top Ni layer while leaving the bottom 
Ni layer free to rotate in a small applied magnetic field. 

Four-probe resistance measurements with the current 
in the plane of the multilayer were performed to deter- 
mine T c . Samples had lateral dimensions 4.3 mm x 1.6 
mm. The T c of each sample was defined to be the temper- 
ature at which the resistance dropped to half its normal 
state value. Fig. 1 shows the results for T c measure- 
ments for samples from several sputtering runs, where d s 
was varied between 16-52 nm. T c shows a strong depen- 
dence on the superconductor thickness close to a critical 
thickness, d / , where the sensitivity to ferromagnetism is 
enhanced. There is no superconductivity above 36 mK 
for d s < df « 16.5 nm. 

The magnetic configuration of our structures was veri- 
fied on simultaneously sputtered samples of larger lateral 
size, in a SQUID magnetometer. Fig. 2 shows a plot of 
magnetization vs. applied field H for a sample with d s = 
18 nm taken at 100 K. The narrow hysteresis loop near H 
— shows the switching behavior of the free Ni layer with 
a coercive field H c — 35 Oe. The wider loop shows switch- 
ing of the pinned layer and is shifted to nonzero H due to 
the exchange bias between the top Ni layer and the FeMn. 
Applied fields of ±100 Oe switch the spin- valve between 
the P and AP configurations. The nearly zero net mag- 
netization observed at -100 Oe indicates very good AP 
alignment between the pinned and free Ni layers, while 
the nearly saturated magnetization observed at +100 Oe 
indicates good P alignment. Similarly good alignment of 
the P and AP states can be achieved at low temperature. 
The inset to Fig. 2 shows a minor hysteresis loop with 
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FIG. 2: Magnetization vs. applied field for a d s — 18 nm 
sample measured at T = 100 K. At ±50 Oe the free bottom 
Ni layer switches while the pinned top Ni layer switches at 
-500 Oe. Inset: minor loop measured at T = 2.29 K showing 
the switching of the free Ni layer. 

H c « 50 Oe taken at 2.29 K, which corresponds to the 
middle of the superconducting transition for this sample. 
We obtain the same behavior for temperatures above and 
below the transition temperature. 

Measurements of T c p and were performed by al- 
ternating the applied field between +100 and -100 Oe, as 
the temperature was slowly decreased through the tran- 
sition region. The largest shift in critical temperature, 
AT C = T ( r p — T c p , should occur in samples with the Nb 
thickness close to d° r . Fig. 3 shows a plot of R vs. T 
for a sample with d s — 17 nm, measured in a dilution 
refrigerator. Two distinct transitions are observed for 
P and AP alignment, with a separation in temperature 
AT C w 28 mK. A second sample with d s = 17 nm showed 
a AT C k 41 mK, but with a slightly broader transition 
centered at 0.34 K. Samples with d s = 18 nm and T c be- 
tween 2 and 3 K exhibit values of AT C of only a few mK, 
similar to the CuNi/Nb/CuNi samples measured previ- 
ously [H [l3|. 

The inset to Fig. 3 shows a plot of R vs. H for the 
first d s = 17 nm sample at a temperature in the middle 
of the transition (0.51 K). The data clearly show well- 
established P and AP states at ±100 Oe, respectively, 
with a difference in resistance of 1.5 £1. Above the tran- 
sition the resistance does not change perceptibly when 
switching from P to AP alignment. An interesting fea- 
ture of the R vs. H curve is the behavior of the resistance 
as the field is swept down from +150 Oe towards -50 Oe 
and as the field is swept up from -150 Oe towards +50 Oe. 
In both cases the resistance increases to a value higher 
than that of the P state after the field passes through 
zero. We believe this behavior involves the breaking of 
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FIG. 3: Resistance vs. temperature for the P and AP states 
of a d a = 17 nm sample measured in ±100 Oe. Two distinct 
transitions are observed, with AT C — 28 mK. Inset: Resis- 
tance vs. applied field at T = 0.51 K (dotted line in main 
graph). 



the free ferromagnetic layer into domains when H ks H c . 
The domain- wall fringe fields penetrate the superconduc- 
tor, thus suppressing T c slightly and producing a higher 
resistance. Note that this effect is opposite to that ob- 
served by other groups 0, H3|, where inhomogeneous 
magnetization led to enhanced superconductivity in F/S 
bilayers. In those experiments, the domain size must 
be smaller than the superconducting coherence length so 
that the Cooper pairs sample multiple domains |2JJ , and 
the magnetic field penetrating into the superconductor 
must be small. 

The critical temperature of F/S/F trilayers in the P 
and AP states has been calculated theoretically by sev- 
eral groups 0, 0, |U HH 13 • Since many experiments 
employ ferromagnetic alloys, the usual approach involves 
solving the Usadel equations in the dirty limit for both 
the superconductor and the ferromagnet. (The dirty 
limit applies to S when Is < £bcs = ^fs7/7r 2 fc_B?co, and 
to F when lp < £ F , where Is and lp are the electron mean 
free paths in S and F, and T c q is the transition tempera- 
ture of the bulk superconductor.) In our case, however, 
the ferromagnetic metal is both pure and strong, thus in 
the clean limit lp > £ F . Hence we use the theory of 01 
as modified in section 3.2 of |25| to make it more appro- 
priate for the clean limit. This theory does not, however, 
incorporate a full description of the majority and mi- 
nority spin bands of a strong ferromagnet, with different 
DOS, vp, and transmission coefficients. The expression 
for the normalized critical temperature of the trilayer is 
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where t c = T c /T c q and T c q is the critical temperature of 
an isolated Nb film of the same thickness as the one in the 
trilayer. The function </> is determined from the condition 
</>tan0 = R for the P state or ((/>tan(/> — R')(R' t&ncf) + 
cf>) — (i?") 2 tan0 = for the AP state, where the complex 
function R = R' + iR" is given by: 



R = 



d s N F v F £c 



1 



& 2N S D S y/l - i£ F /l F + 2/T F 
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Eq. (2) is valid when the ferromagnets are thick enough 
so that the tanh functions in [25| can be set to 1. This 
assumption is validated by data on Nb/Ni bilayers 0] 
where oscillations in T c (d F ) are completely damped for 
d F > 4 nm. The dimensionless parameters that en- 
ter into this theory are the ratios d s /£s, £,f/If, the 
S/F interface transparency T F , and the combination 
N F v F £s/2NsDs- N F and Ns are the densities of states 
at the Fermi energy of the F and S layers, v F is the Fermi 
velocity of the ferromagnet, and Ds is the diffusion con- 
stant of the superconductor. 

To avoid fitting the data with four free parameters, we 
follow the strategy outlined by Lazar et al. § and by 
Sidorenko et al. |l6| . We determine the superconducting 
coherence length, £g, from measurements of the critical 
field vs. temperature of isolated Nb films, with the mag- 
netic field applied perpendicular to the film plane. For 
films in the thickness range 20-50 nm, the values of £s 
are close to 6 nm, which we use for our fits [2(|. From 
the asymptotic form of Eq. (1) as t c — > 0, one finds 
20 2 /(C/Cs) 2 = 1/47, where 7 = 1.781. Substituting 
df 16.5 nm and using Eq. (2) (while ignoring the 
small imaginary term), we obtain the constraint 



N F v F £ s 



1 



r tan0 c 



0.24 (3) 



(1) 



2N s D s (d?) 1 + 2/T F (dy/&) 

Estimates of the product N F v F vary substantially in the 
literature. From and 0], we obtain respectively 
N F = 1.77 • 10 48 J-'m" 3 and v F = 0.28 • 10 6 m/s. 
Fierz et al. [2^, however, quote p F l F — 0.7 — 2.3 ff2m 2 
for Ni, which when combined with the Einstein relation 
l/p F l F = N F v F e 2 /3 imply values 3-10 times smaller for 
N F v F . Combining these values with Ns = 5.31 • 10 47 
J _1 m -3 [2jj and using our measured Ds(df) — 2.8 • 
10~ 4 m 2 /s, we obtain T F = 0.05 - 0.6. The bulk resis- 
tivity of our sputtered Ni films at 4.2 K is p F = 33 nfim, 
which leads to values of I F between 7 and 70 nm, given 
the range in p F l F quoted above. Since the Ni used in our 
trilayers is thin, lp is probably limited by surface scat- 
tering, so we use the lower estimate lp = 7 nm, hence 
£p/l F « 0.1. In fact, the fit to T c (d s ) is quite insensitive 
to the values of Tp and £, F /l F . We used £ F /lp = 0.1 and 
T F = 0.3 to obtain the curve shown in Fig. 1, which fits 
the data remarkably well. 

A more stringent test of the theory is the prediction of 
AT C , which depends sensitively on both Tp and S, F /l F . 
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FIG. 4: Symbols: AT C vs. T c for our 11 thinnest samples. 
The line represents a fit using £f/If = 0.7 and Tf = 1-0, 
values larger than our best estimates. 



Thickness deviations from nominal values produce scat- 
ter in plots of T c or AT C vs. d s , therefore Fig. 4 shows a 
plot of AT C versus T c . If we calculate AT C using our best 
estimate of £f/If and the upper limit of Tp given above, 
the maximum value of AT C is only a few mK when T c 
is well below IK- hardly visible on Fig. 4. If we relax 
the constraints we have placed on the parameters, and 
instead try to produce the best fit to the AT c (d s ) data, 
we find that a reasonable fit can be obtained when £p /If 
is allowed to be much larger than our original estimate. 
Fig. 4 shows a fit using £f/If — 0.7 and Tp = 1.0. Sim- 
ilar curves can be produced by simultaneously varying 
£,f/If and Tp while keeping their product nearly con- 
stant. Fitting the AT C data requires letting £,f/If exceed 
our estimate substantially. Our lp estimate may be too 
large, because the resistivity is dominated by the longer 
of the majority or minority band If, whereas the F/S 
proximity effect depends on the shorter of the two |8j ■ A 
shorter / p is also implied by the observation of complete 
damping of T c oscillations in Nb/Ni bilayers for dp > 4 
nm [16j. Nevertheless, producing a reasonable fit to our 
AT C data entails either increasing £,f/If beyond the clean 
limit, or increasing Tp beyond our original estimate. 

In conclusion, we have observed a large difference in 
T c between the P and AP magnetic states of Ni/Nb/Ni 
trilayers, with T c p < T^ F '. Recently, Ruzanov et al. 
pjof reported a T c difference between the P and AP 
states of a Ni .sFeo.2/Nb/Ni .8Feo.2 trilayer, but with 
T P > T^ p . Understanding these opposing behaviors in 
F/S systems with strong ferromagnets will require fur- 
ther experiments, as well as theoretical models able to 
account for the complexity of real ferromagnets [3l|. 
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